Search for neutrinos from Gamma-Ray Bursts with 

ANTARES 



Eleonora Presani 

Nikhef, Science Park 105, 1098 XG Amsterdam 

Abstract. A method to search for neutrino induced showers from gamma-ray bursts in the ANTARES detector is presented. 
ANTARES consists of a three-dimensional array of photosensitive devices that measure Cherenkov light induced by charged 
particles produced by high energy neutrinos interacting in the detector vicinity. The shower channel is complementary to the 
more commonly used upgoing muon channel. The corresponding detection volume is smaller, but has the advantage of being 
sensitive to neutrinos of any flavour. 
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INTRODUCTION 

Through this analysis, the Antares experiment is sensitive to the electromagnetic and hadronic showers induced in the 
detector by the interaction of neutrinos originating from Gamma-Ray Bursts. 

In the Fireball model, neutrinos are produced during the GRB explosion via photo-pion production in collision of 
ultrarelativistic protons with photons in the jet of the GRB. Although the production mechanism favours the muon 
neutrino it can safely be assumed, due to neutrino oscillation, that the flavour ratio at Earth is v e : : V T = 1 : 1 : 1 . 

While the standard track search is sensitive only to muon neutrinos, interacting via charged current interaction 
with nuclei, a shower analysis is sensitive to all flavours. In order to deal with the high background introduced by 
atmospheric muons, a track search is generally limited to upgoing neutrinos, i.e. neutrinos that traversed the Earth and 
interacted near the detector. The current shower analysis takes advantage of the time coincidence with a triggered GRB 
to reduce the background and can therefore extend the field of view to the entire solid angle. To estimate the discovery 
potential and the sensitivity of the detector, a Waxman-Bahcall flux, averaged for a single GRB, was assumed. The 
analysis consists of the search for coincidences between well reconstructed showers in Antares and photons triggered 
in one of the satellites which subscribed to the IPN (Interplanetary Gamma-Ray Burst Timing Network). 



NEUTRINO PRODUCTION IN GRBS 

In the Fireball model [1] charged particles can be accelerated to high energies via the Fermi acceleration mechanism. 
Neutrinos are primarily produced via the p + y (photo-hadronic) process. The Waxman-Bahcall flux for a single 
average Gamma-Ray Burst [2] was used for the shape and normalization of the flux distribution. It is a diffuse flux, and 
it is therefore necessary to derive the flux of a single average GRB. This procedure relies on a number of assumptions: 
the average duration for each GRB is given by T90=50 s and 667 GRBs are, on average, detected per year by satellites. 
This number is based on the observations of the CGRO mission. The resulting flux can be written as: 
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where the 4n takes into account the solid angle in which the flux is generated, N y J" r is the number of seconds in a year, 
and the denominator represents the average number of GRB multiplied by their average duration (chosen to be 50 s). 



THE ANTARES DETECTOR AND EVENT SIGNATURES 



The Antares (Astronomy with a Neutrino Telescope and Abyss environmental RESearch) neutrino telescope [3] is 
located at around 40 km off Toulon, at a depth of 2475 m. The photosensors are arranged along a length of <~ 350m 
of vertical cable. A total of 885 optical sensors are carried by 12 flexible lines. Each line comprises up to 25 detection 
storeys, each equipped with three downward-looking 10-inch photomultipliers (PMTs), oriented at 45° relative to the 
vertical. Each PMT is housed in an Optical Module (OM) that consists of a 17-inch glass sphere in which the optical 
connection between the PMT and the glass is ensured by an optical gel. 

The majority of neutrinos that reach the Earth will just pass through it. However, it is possible for them to undergo 
a weak interaction with a nucleon. Due to the extremely small cross section, a very large target mass is necessary to 
attempt the detection. The reaction that is used in high energy neutrino detectors to reveal this particle is the neutrino 
deep inelastic scattering (DIS) with a matter nucleon. If a W boson is exchanged (charged current interaction), a 
charged lepton and an hadronic shower are generated, while if a Z° boson is exchanged (neutral current) an hadronic 
shower is followed by another neutrino. When a neutrino interacts in the vicinity of the detector via a charged current 




FIGURE 1. Schematic view of a "shower" like event in the Antares detector. Neutrinos of any flavour will produce an hadronic 
shower in the detector when interacting via NC (left). Electron neutrinos will produce a shower in the detector for both NC and CC 
interactions (right). 

(CC) interaction, it generates a relativistic charged particle that will travel through the sea water in the detector. When 
the speed of a charged particle exceeds the velocity of light in the medium through which it is travelling it will induce 
the emission of a radiation called Cherenkov radiation [4]. The light emission creates a wavefront where the emitted 
light is coherent. The wavefront forms a cone with its apex at the travelling particle. The opening half-angle of the 
cone C can be written as cos C = 1/fin. The sea water at the location of the ANTARES neutrino telescope has a value 
of n of about 1.35, thus the value of the Cherenkov angle is about 42.2°. 

Neutrinos of any flavour interacting via a NC will produce a neutrino, which remains undetectable, and an hadronic 
shower, consisting of many interaction products. Due to the small volume where the shower takes place the shower 
appears in the detector as a point-like light source developing isotropically in time. Tau neutrinos can also be detected 
with shower analysis, as their interaction and decay produce showers. A neutrino that undergoes a NC interaction will 
not produce a charged lepton capable of leaving a detectable track (Fig. 1). If the main vertex of the interaction is 
in the vicinity of the detector, it is possible to see light as if being emitted from a point source. For an hadronic (or 
electromagnetic) shower to be measured, it is necessary for the primary vertex to occur in or around the instrumented 
detector volume. This reduces the probability of production of a reconstructible shower. 



The majority of the photons measured in Antares (called "hits") are due to optical background and to Cherenkov light 
induced by down-going muons created in interactions of high energy cosmic rays in the atmosphere. For this work a 
strict hit selection is used to distinguish between background and signal hits. The hit selection starts by searching for 
clusters of hits in time, because these are more likely to be due to a signal than background. During the hit selection, 
the detailed geometry of a single storey is ignored and all hits on a floor are considered together to form clusters. For 
each event, all hits on one storey are time ordered, all hits that occur within 20 ns of each other are merged together. 
Using these as seeds, all the hits that are causally connected to the clusters are kept. To be sure that only signal hits are 
selected, only hits compatible with being generated by the same particle are selected [5]. Given a shower occurring 
at a position (X true ,Y true ,Z true ) and time t true and emitting light, the expected arrival time, t' exp , of photons is given by 
Pythagoras' rule: 



SHOWER EVENT RECONSTRUCTION 




(2) 




FIGURE 2. Time and spatial resolution of the reconstruction. The first shows the difference between the reconstructed and the 
true time of the shower. The time resolution of the reconstruction taken to be the RMS of the distribution is 5.1 ns. The second 
shows the logarithm of the three dimensional distance between the true and reconstructed shower vertex. The spacial resolution is 
best for electron neutrinos, with a median value of 3.0 m. 
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FIGURE 3. Reconstruction efficiency as a function of neutrino energy for an average GRB flux. See text for details. 



where n is the group refractive index of sea water (n = 1.3797) and c is the speed of light in vacuum. To reconstruct 
the time t reco and position (X reco ,Y reco ,Z reco ) of the shower, these unknown parameters are varied until the difference 
between the expected time and the measured time of each selected hit is minimised. The minimising function chosen 
for this analysis is the so called M-estimator function, defined as: 



The errors a, could in principle be different for each hit. In the Antares experiment, all optical modules are observed 
to have a similar resolution and, therefore, the same value of a, ■ = o = 1 ns is used for each hit. 

The resolution of the reconstruction is shown in Figure 2, where spectrum is generated weighting the events with 
a Waxman-Bachall spectrum.. The time resolution is measured to be 5.1 ns. The right plot in Figure 2 shows the 
three dimensional distance between the true and reconstructed shower. The spatial resolution of the reconstruction is 
given by the median of this distribution. It is better for electron neutrinos, with a value of 3.0 m, and worse for muon 
neutrinos, with a value of 9.3 m. The reason for this difference is that electron neutrinos always produce a shower-like 
signal. On the other hand, only muon neutrinos that undergo a NC interaction will produce a pure shower signature. 
The second peak, especially visible for muon neutrino induced events, is due to Bremsstrahlung and electromagnetic 
showers generated along the muon track that are reconstructed far away from the interaction vertex. The efficiency 
of the reconstruction has also been studied. The efficiency of the shower reconstruction itself depends on the criteria 
used to select reconstructed showers. A minimum bias reconstruction efficiency can be determined by accepting all 
events in which the reconstruction successfully found a minimum. For this work, in addition, 5 selected hits on at least 
2 lines were required for the minumum selection. The efficiency of the reconstruction as a function of the neutrino 
energy after this first selection is shown in the first plot of Figure 3. The reconstruction is considerably less efficient 
at low energies, where showers produce fainter light and less hits. At high energy practically all triggered events are 
reconstructed. The right plot in Figure 3 shows the efficiency with which the shower fit procedure determines very 
accurately the position and time of the shower. Only events for which the reconstructed shower was within 3 m and 
10 ns of the true shower vertex are used for calculating this efficiency. Even with these strict cuts on the quality of 



the reconstruction, the efficiency stays between 30% and 50%, especially in the energy range on which the analysis is 
focused. 



NEUTRINO INDUCED SHOWER ANALYSIS 

The concept of time correlated analysis is very simple: given a certain search time window, one looks for a neutrino 
event in the same time window as the one of a gamma-ray burst event. Quality cuts on the reconstructed events 
are tuned in order to reduce the background and to optimize the discovery potential. Since this analysis is focused 
on shower events, no angular cut is applied, and the entire sky is considered. The main advantage of using a time 
correlated search is the efficient background rejection, due to the small time window considered. The time window 
used for the quality cut optimization is 100 seconds. Each GRB will be analysed for the duration of its T90. A set of 
data corresponding to 295.8 days of data taking during the year 2008 has been used to extract the background level. 
During the same year 65 GRBs have been recorded and considered for this work. 

The quality cut optimization was done on three parameters: the quality parameter (the M-estimator of the recon- 
struction), the number of hits assumed to come from the shower and the number of lines used for the reconstruction 
of the shower. It is important to have well-reconstructed shower, therefore the quality parameter is very significant. 
The hits coming directly form the shower are the set of all the selected hits that have a time residual smaller than 
15 ns. This parameter is important as it gives some information on how good the shower hypothesis is for that event. 
If the event does not behave as a shower, the hits will not be distributed as expected. The number of detector lines 
used in the reconstruction is correlated with the energy of the shower. The optimization of the quality cuts is made 
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FIGURE 4. Sensitivity of the Antares detector for neutrino induced showers from GRBs. 

in order to maximize the discovery potential of the analysis. The optimized values lead to a 5(7 discovery if 4 events 
are measured within the time window of the GRB (100 s in the optimization procedure) or a 3(7 discovery when 3 
events are measured. This method will then be applied for each observed GRB. The 90% C.L. average upper limit of 
Antares (or sensitivity) for a neutrino induced shower in coincidence of GRB trigger is shown in Figure 4. The dashed 
line shows the model used for the estimate of the neutrino rate [2], The thick black line is the average upper limit 
on the neutrino flux for 2008 data, calculated using the Feldman-Cousins method [6]. The sensitivity of Antares for 
showers in coincidence with a GRB is Ey ^ < 1.67 x 10~ 2 GeVcm- 2 s _1 for 10 s < E v < 10 7 . This is the sensitivity 
corresponding to the 2008 data set, during which the average background rate was 1.782 x 10~ 4 Hz. 
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